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THE  ATTACHED  FLOW  OVER  A  SLENDER  DELTA  WING  WITH  LEADING-EDGE  FLAPS 

by 

A.  J.  Peace* 


SUMMARY 

The  flow  over  a  delta  wing  with  leading-edge  flaps  is  considered 
within  the  framework  of  slehder-body  theory.  To  avoid  flow  separation  at 
the  leading  edges  of  the  flaps  a  Kutta  condition  is  applied  there  which 
uniquely  determines  the  angle  of  attack  for  a  given  flap  configuration. 
The  lift  coefficient  and  lift-dependent  drag  factor  are  then  calculated 
for  this  condition  and  compared  with  the  prediction  of  slender  thin-wing 
theory,  in  which  the  wing  boundary  condition  is  fully  linearized. 


*  School  of  Mathematics  and  Physics,  University  of  East  Anglia.  This  work 
was  undertaken  during  a  period  at  RAE  under  the  CASE  scheme. 
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1  INTRODUCTION 

The  idea  of  deflecting  downwards  the  outboard  part  of  a  delta  wing  about  a  hinge- 
line  or  shoulder-line  near  the  leading  edge  has  been  current  for  a  considerable  period  *  . 
It  has  been  studied  both  as  a  simple  form  of  built-in  conical  camber,  intended  to  recover 
some  part  of  the  axial  force  which  is  lost  through  separation  at  the  leading  edge  of  a 
plane  slender  wing;  and  as  a  variable-geometry  device,  intended  to  enlarge  the  range  of 
lift  coefficient  over  which  an  acceptable  performance  can  be  achieved.  Many  measurements 
have  confirmed  that  lower  drag  at  lifting  conditions  can  be  obtained  with  such  a  leading- 
edge  flap. 

The  benefit  arises  because  the  pressure  differential  across  the  wing  is  largest 
close  to  the  leading  edge,  so  that  a  downward  deflection  of  the  surface  produces  a  force 
in  the  direction  of  flight.  Since  the  downward  deflection  also  reduces  the  size  of  the 
peak  in  the  pressure  differential,  guidance  is  needed  in  the  choice  of  an  appropriate 
form  and  amount  of  deflection  for  a  particular  flight  condition.  Providing  this  guidance 
is  not  easy  because,  in  general,  the  flow  will  separate  at  the  leading  edge  of  the  wing 
to  form  a  vortex. 

However,  if  the  angle  of  incidence  is  not  too  large,  it  is  possible  to  droop  each 
element  of  the  leading  edge  in  such  a  way  that  it  points  directly  into  the  local  on¬ 
coming  flow,  so  that  any  tendency  of  the  flow  to  separate  at  the  leading  edge  is 
suppressed  at  this  condition.  In  particular,  if  the  flow  past  a  delta  wing  is  conical, 
this  condition  can  be  brought  about  by  deflecting  a  leading-edge  flap  downwards  about  a 
straight  hinge-line,  the  amount  of  deflection  being  greater  for  larger  angles  of  incidence. 
For  a  wing  at  this  'attachment  incidence'  the  flow  can  be  calculated  by  the  fully- 
linearized  theory  of  attached  flow,  and  it  emerges  that  values  of  the  lift-dependent  drag 
not  far  above  the  minimum  (corresponding  to  elliptic  loading)  are  predicted  for  suitably 

designed  wings.  It  also  emerges  in  particular  cases  from  both  experiment  and  calcula- 
4  . 

tion  that  it  is  beneficial  to  use  rather  less  leading-edge  droop  than  would  be  required 
to  achieve  the  condition  of  attached  flow. 

Despite  the  developments  that  have  taken  place  in  modelling  the  separated  flow, 
based  on  slender-body  theory  and  on  panel  methods,  there  is  no  adequate  method  for  cal¬ 
culating  the  lift  and  drag  of  wings  with  small  leading-edge  vortices  at  subsonic  speeds. 

It  is  therefore  useful  just  to  know  how  much  droop  of  a  particular  form  is  needed  to 

suppress  the  leading-edge  separation.  For  the  modest  amount  of  built-in  camber  appro- 

3  5 

priate  to  a  cruising  condition,  it  appears  ’  that  a  fully-linearized  theory  predicts 
this  adequately.  However,  if  a  leading-edge  flap,  capable  of  deflection  through  large 
angles,  is  provided  to  improve  performance  at  the  much  larger  values  of  the  lift 
coefficient  associated  with  take-off  and  manoeuvre,  it  is  not  clear  that  the  fully- 
linearized  theory  will  be  adequate,  particularly  as  calculations6  for  a  smoothly  cambered 
wing  revealed  differences  between  the  predictions  of  slender-body  theory  and  the  (fully- 
linearized)  slender  thin-wing  theory.  The  present  calculations  were  undertaken  to  dis¬ 
cover  how  much  deflection  of  a  leading-edge  flap  is  required  to  suppress  leading-edge 
separation  on  a  flat  slender  delta  wing,  and  what  lift  coefficient  and  lift-dependent 
drag  factor  are  predicted  at  this  condition,  and  to  assess  how  well  these  are  predicted 
by  the  fully-linearized  theory. 
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In  the  slender-body  approximation  the  velocity  potential  satisfies  Laplace's  equa¬ 
tion  in  two  dimensions  in  each  transverse  or  'cross-flow'  plane.  The  solution  of  this 
equation  is  achieved  using  complex  variable  techniques.  He  regard  the  cross-flow  plane 
as  being  a  conplex  plane  and  use  a  conformal  transformation  to  map  the  region  of  the 
right-hand  side  of  this  plane  outside  the  cross-section  of  the  wing  into  the  upper  half 
of  the  transformed  plane  so  that  the  half-wing  is  mapped  onto  the  real  axis  of  the  trans¬ 
formed  plane. 

The  geometrical  parameters  defining  the  wing,  ie  the  flap  length  and  deflection 
angle,  lead  to  the  aerodynamic  properties  of  the  configuration.  The  results  are  obtained 
by  prescribing  these  two  parameters  and  then  calculating  the  required  angle  of  attack  to 
produce  an  attachment  line  on  the  leading  edge:  the  lift  and  drag  then  follow. 

2  STATEMENT  OF  THE  PROBLEM 

A  cartesian  coordinate  system  with  origin  at  the  apex  of  the  wing  is  introduced. 

The  x-axis  lies  along  the  wing  centre  line  with  the  y-axis  running  to  starboard  and  e-axis 
vertically  upwards,  so  that  the  set  Oxyz  is  right-handed.  The  local  semi-span  is  given 
by  s  »  Kx  ,  where  K  ■  tan  y  ,  and  2y  is  the  apex  angle  of  the  wing  with  flap 
undeflected.  With  a  flap  deflection  0  in  the  cross-flow  plane  x  -  constant  ,  the  con¬ 
figuration  of  Fig  1  is  obtained.  The  apex  angle  of  the  undeflected  part  of  the  wing  is 
2y'»  with  the  flap  shoulder  a  distance  ns  from  the  wing  centre  line;  n  is  thus  a 
conical  coordinate  of  the  flap  hinge-line.  The  flap  length  is  therefore  (1  -n)s  and 
nK  ■  tan  y'  .  The  wing  is  shown  in  Fig  1  at  incidence  a  to  the  undisturbed  stream  of 
velocity  U  ,  which  is  resolved  into  components  U  cos  a  along  Ox  and  U  sin  a  along 
Oz  .  These  components  are  approximated  by  U  and  Ua  respectively.  Hence,  we  seek  a 
function  $  such  that  the  velocity  potential  4  is  given  by 


®  ■  Ux  ♦  4  • 


0) 


Note  that  the  extra  term  in  the  slender-body  potential  which  takes  into  account  the  thick¬ 
ness  of  the  configuration  vanishes  here  as  we  are  considering  wings  with  no  thickness. 


Slender-body  theory  shows  that  the  function  4 
Laplace  equation 

2  2 

+  .  o 

2  2 

3y  3  s 


satisfies  the  two-dimensional 


in  each  cross-flow  plane.  We  are  therefore  treating  the  transverse  velocity  field  as 
though  the  flow  were  two-dimensional,  the  three-dimensional  nature  of  the  problem  only 
entering  through  the  boundary  conditions,  which  are  that  the  velocity  normal  to  the  wing 
is  zero,  ie 

n  •  V*  -  0  (2) 


where  n  is  the  unit  outward  normal  to  the  wing,  and  that 


$  -  Uza 


at  large  distances  from  the  wing. 
He  introduce  new  variables 


A  A  A 

and  note  that  $  satisfies  Laplace's  equation  in  the  (y,z)  plane.  He  therefore  seek  a 
complex  potential 

H(Z)  -  $(Z)  +  i$(Z)  (5) 

as  ^  ^ 

where  tp  is  the  stream  function  corresponding  to  the  potential  $  ,  and  Z  -  y+iz  is  a 
non-dimensional  complex  coordinate  in  the  cross-flow  plane.  The  solution  is  obtained  by 
conformally  mapping  the  right-hand  half  of  the  Z-plane  exterior  to  the  cross-section  of 
the  wing  onto  the  upper  half  of  the  t-plane  so  that  the  contour  A^BCDEFA^  is  mapped 
into  the  real  axis  of  the  t-plane  (Fig  2),  The  potential  $  can  then  be  found  in  the 
t-plane  as  that  due  to  a  distribution  of  sources  on  the  real  axis,  together  with  a  uniform 
flow,  and  inversion  of  the  transformation  will  give  the  potential  in  the  Z-plane. 

3  THE  TRANSFORMATION 

He  denote  the  image  of  L  in  the  Z-plane  by  L  in  the  t-plane.  He  may  specify 
that  t  »  0  at  5  and  denote  the  values  of  t  at  5,  C,  E,  F  by  b,  c,  e,  f  where 

b<c<0<e<f.  (6) 

The  required  transformation  is  of  the  Schwartz-Christof fel  type,  whose  general  form 

||  -  An(t-t.)ni  , 

where  t.  are  the  images  of  the  vertices  of  the  polygon  in  the  Z-plane  and  the  n.  are 
related  to  changes  in  direction  of  the  contour  at  the  vertices.  Inserting  the  appropriate 
values  of  m  and  t.  gives 


dZ 

dt  "  A 


=====(— )6 

'(t  -b)(t  -f)  \t  '*/ 


where  A  is  a  constant.  Then  by  substituting  particular  values  of  t  and  Z  we  obtain 
the  following  equalities: 


n  -  a  F(t)dt  , 


0  -ri)e'lS  -  A  f  F(t)dt  , 
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(1  -n)e-i<e+,r) 


e 

A  J  F(t)dt  , 
0 


(10) 


-  n 


A  J  F(t)dt  . 
e 


(11) 


In  evaluating  the  right-hand  sides  of  equations  (8)  to  (11)  we  nust  ensure  the 
appropriate  value  is  taken  in  the  detenoi nation  of  the  integrands  as  they  are  not  single— 
valued  functions.  Writing  t  ■  Itle1®,  t-b  ■  |t-b|e  **  etc,  and  A  ■  |A|e  a  ,  we 
define 


dZ 

dt 


i(e.8,-ieb-lvlvK) 


Then  (8)  gives: 


A  f  - . - (l—f*  dt 

1  l  /(t  -b)(f  -t)  '*  "  ' 


(12) 


As  the  left-hand  side  of  (12)  is  real  and  positive  and  the  integral  is  real  and  positive 
we  oust  choose  A  to  be  imaginary  and  negative  to  give  the  correct  rotation  between  the 
Z-  and  t-p lanes.  It  is  possible  to  choose  | A |  -  1  ,  so  preserving  the  boundary  condition 
at  infinity,  because  there  is  still  some  freedom  in  the  transformation. 


Hence 


and  therefore 


A  •  -  i 


c 

/ 


- 1 

(1ZS. 

/(t-b)(f-t) 

\e-t, 

of  the  radicals  are  : 

rt-cX1 

/(t-b)(f-t)  1 

<*-t) 

6  | 

't  -cl1 

/(t-b)(f -t) 

ke  -t / 

c  ( 

't  -c\ 

/(t-b)(f  -t) 

kt  -  a/ 

B/n 


dt  ■  n 


B/n 


dt  -  1  -  n  , 


B/n 


dt 


1  -  n  i 


B/n 


dt 


(13) 


(14) 


(15) 


(16) 


(17) 


We  now  have  in  (14)  Co  (17)  four  equations  in  four  unknowns,  b,  c,  e,  f,  depending 
on  the  two  parameters  n  and  0  ,  giving  a  problem  which  is  apparently  determinate  and 
capable  of  numerical  solution. 

However,  before  numerical  treatment  is  possible  the  singularities  in  the  integrands 
must  be  removed.  To  perform  this  we  use  the  following  substitutions: 


in  (14)  put 

t  -  b 

- 

2 

r  , 

in  (16)  put 

e  -  t 

s 

k 

r  , 

where  k  =  1/(1  -  (B/ir)), 

and  in  (17)  first  put 

f  -  t 

= 

2 

r  , 

and  then  put 

f  -  e  -  r2 

= 

_k 

r  » 

for  i  <  r  <  (f  -  e)  ^ , 

■  interval  has  to  be  split. 

Equations  (14)  to  (17)  can 

be  written 

in 

the  form: 

F(x) 

0  , 

(18) 

where  F  is  a  vector  of  four  nonlinear  functions  of  the  vector  x  *  (b,c,e,f)  .  This 
system  is  solved  by  a  multi-dimensional  Newton  scheme.  The  solution  of  (18)  for  the 
complete  range  of  ri  and  8  is  most  easily  obtained  by  keeping  ii  constant  and  spanning 
the  range  of  B  .  The  solution  for  a  pair  (n,B)  can  then  be  used  as  an  initial  guess  at 
the  solution  for  a  neighbouring  pair.  The  region  of  the  (n,B)  parameter  space  which  is 
physically  and  numerically  accessible  is  discussed  in  section  5. 

4  SOLUTION  OF  THE  PROBLEM 

4. 1  Boundary  conditions 

A 

The  boundary  condition  on  the  undeflected  part  of  the  wing  requires  that  ■*  0 
on  that  part  of  the  contour.  From  the  geometry  in  Fig  1  it  can  be  shown  that  the  unit 
outward  normal  to  the  flap  is 

n  -  -  g  ,  -  sin  6  sin  y',  sin  6  cos  y',  cos  8  cos  y'J  .  (19) 

/l  -  cos^B  sin^y* 


We  require  that  the  velocity  normal  to  the  flap  is  zero,  so,  by  (1),  (2)  and  (19),  the 
velocity  normal  to  the  contour  in  the  cross-flow  plane  is 


an 


+  U  sin  8  sin  y' 

/l  -  cos  ^6  sinZY' 


(20) 


where  the  ♦  and  -  refer  to  the  upper  and  lower  surfaces  of  the  flap,  >spectively. 


and  n  is  the  distance  normal  to  the  contour 


8 


Therefore  Che  required  non-dimensional  normal  velocity  in  the  t -plane  is 


l£ 

.  iil 

1  — 1 

-  lii |<£l 

* 

3n  t 

35  lz 

1  dt  J 

U  3n  J  dt  I 

where  n  is  the  non-dimensional  distance  normal  to  the  contour,  and  higher-order  terms 
in  K  have  been  neglected  in  the  expansion  of  the  right-hand  side  of  (20),  and  here  the 
+  and  -  refer  to  the  line  segments  CD  and  DE  ,  respectively.  The  normal  velocity  is 
zero  elsewhere  on  the  contour. 

To  realise  this  normal  velocity  we  introduce  a  continuous  non-dimensional  source 
distribution  in  the  t-plane  with  strength  per  unit  length  at  t'  given  by 


c  <  t’  <  0, 

0  <  t'  <  e,  ►  (22) 

elsewhere. 


2nK  sin  e  yr  , 
dt  f 


-  2nK  sin  6  , 


We  must  also  satisfy  the  boundary  condition  (3). 
4.2  Complex  potential 


As  stated  in  section  2  we  shall  form  a  complex  potential,  W  ,  to  Batisfy  the 
boundary  conditions  of  section  4.1. 


Wj  ,  where 


For  the  contribution  from  the  source  distribution  (22)  we  obtain  a  complex  potential 
'here 

dW,  f  i  i  j.  i  « 


■df  ■  2nK  sin  0  j  ±  ±  |g|  T±rr  dt* 


where  <S f  is  the  curve  CE  in  the  t-plane,  and  the  +  and  -  refer  to  the  portions  of 
the  curve  stated  previously. 


To  satisfy  (3)  we  require 


So  (24)  and  (13)  imply  that 


dW  ~  _ 
dZ 


dW  _  _ 
dt 


as  | Z )  ®.  (24) 


as  1 1 1  -*• 


^  -  a  +  f(t)  (25 

where  |f(t)|  -*■  0  as  |t|  -*■  «  .  The  confined  potential  which  satisfies  (25)  and  gives 
the  correct  normal  velocity  as  produced  by  (23)  is 
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or 


dW 

dt  “ 


/  0 

fjK  Bin  6  I  f 

’  \r 


i 


dt' 


U 


at  ♦  SUlHLl 


in(t-t')dt'  . 


(26) 


(27) 


4.3  Attachment  condition 


To  avoid  separation  at  the  leading  edge  o£  the  flap  we  impose  a  Kutta  condition  to 
remove  the  singularity  in  velocity  and  pressure.  This  will  determine  the  value  of  a 
(afl,  say)  that  gives  attached  flow  for  a  given  flap  configuration.  For  attachment  we 
require  that  the  velocity  be  finite  at  the  leading  edge,  which,  since  dZ/dt  =  0  there, 
leads  to 


dW 

dt  t-0 


-  0 


(28) 


Therefore  from  (26)  and  (7)  we  obtain  (writing  t'  ■  r) 


r?K  sin  6 
ir 


o  e\ 

Vc  0 


/(? 


==(—) 

-b)  (r  -  f)  'r~e' 


B/ti 


-r 


dr  , 


and  simplifying. 


a  n  sm 
K 


in  B  f  _ I _  fr  -c\ 

"  J  /(r-b)(f-r)  'e_r' 


r_-c\B/lT  dr  , 


(29) 


(30) 


where  again  the  singularity  in  the  integrand  must  be  removed  by  substituting  e-r  -  r 
before  numerical  treatment  can  be  applied. 

We  may  also  note  at  this  point  that  a  weaker  singularity  in  the  velocity  and 
pressure  exists  at  the  shoulder  of  the  flap  and  for  moderate  or  large  flap  deflections 
we  should  expect  separation  at  this  point.  However,  this  feature  is  being  ignored  in  the 
present  analysis  and  the  potential  flow  calculated  here  predicts  smooth  flow  around  the 
shoulder  with  infinite  velocity, 

4.4  Lift 

The  lift  can  be  obtained  conveniently  by  the  use  of  a  distant  control  surface.  If 
the  non-dimensional  complex  potential  W  is  expanded  for  large  values  of  itd  non- 
dimensional  argument  Z  in  a  plane  of  constant  x  ,  we  obtain  an  expression  of  the  form: 


W  -  a,Z  ♦  aQ  +  a_,z"1  +  ...  .  (31) 

The  lift,  L  ,  acting  on  the  portion  of  a  thin  wing  upstream  of  the  plane  of  constant  x 
is  then  given  in  terms  of  the  oefficier  a_j  by  a  particular  case  of  equation  (9.7.11) 
of  Ref  7  as 


10 


(32) 


L  -  2wpU282>{a_I}  . 

Expansion  of  (7)  for  large  t  ,  followed  by  integration,  shows  that 
Z  -  -  i|t  +  Yjj  ♦  Y  In  t  +  Y-1C  *  +  ?  } 


(33) 


.  -  b  eg  eg  f 

where  y  -  -f'  T~2  * 


By  a  similar  argument  to  Andrews8  in  his  Appendix  A  it  can  be  shown  that  the  coefficient 
of  in  t  in  the  expansion  of  Z  (or  that  of  1/t  in  the  expansion  of  dZ/dt  )  must 
vanish.  Hence, 


b.cg  +  eg  f  .  0  # 

2  IT  IT  2 


(34) 


Using  (34)  we  have 


J-f  (e-c)2  -  -fa  (e  -c)  (e  +  c  -  b  - f)  -  j  (b-f)1 

2tt 


(35) 


To  invert  (33)  we  write 


l|cjZ  +  Cq  +  c_|Z  +  0(Z  )| 


whence 


C1  "  1  •  C0 


-r 

*  c-i  '  Y-i 


(36) 

(37) 


Expanding  (27)  for  large  t  gives 


0  e\ 


W  -  -  at  + 


nK  sin  B 


J-j  mt.<r • 


(38) 


but  from  (15)  and  (16) 


0  e\ 


/-  /  M-"'  '  0 


(39) 


so  introducing  (36)  and  (37)  into  (38)  and  using  (39)  gives 


0  e\ 


*•(* ♦  it. * *.,o * if- f  * o<r2) 


(40) 
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Hence 


(41) 


and  so,  at  the  attachment  incidence,  a  , 

’  a 


* 

1  0  f\ 

2  2 

L  -  2irs  pU  < 
a 

.  nK  sin  0 

-  a  y  .  +  - 

a  -1  it 

(42) 


.  .  2 
A  lift  coefficient,  ,  based  on  the  kinetic  pressure,  jpU  ,  and  the  planform  area 

a  2 

of  the  wing  with  flap  undeflected,  s  /K  ,  is  given  by 


(43) 


or  simplifying  the  integral  (writing  t*  =  r) 


4ir 


V- 


1  rj  sin  S 


/ 


v'(r-b)  (f  -  r) 


(44) 


where  y_j  is  given  by  (35).  Here  also  the  singularity  in  the  integrand  must  be  removed 
for  numerical  purposes,  as  before. 

In  the  Appendix,  equation  (44)  is  expanded  for  small  values  of  0  ,  and  it  is 
verified  that  the  leading  term  agrees  with  the  expression  obtained  from  the  fully- 
linearized  slender  thin-wing  theory,  in  which  the  wing  and  flap  boundary  condition  is 
applied  in  the  plane  z  »  0  . 

4.5  Drag 

The  drag  is  obtained  by  first  calculating  the  normal  force  on  the  flap.  The  overall 
drag  can  then  be  calculated  by  resolving  the  normal  forces  on  the  wing  and  flap  as  shown 
later.  The  normal  force  on  the  flap  is  evaluated  by  integrating  the  pressure  jump  across 
it:  the  pressure  is  obtained  from  the  velocity  components  on  the  flap.  The  entire  force 
comes  from  the  pressure  integral  as  there  is  no  localised  suction  force  on  the  leading 
edge  as  we  are  considering  the  attachment  incidence  only. 

If  u  and  v  are  the  non-dimensional  velocity  components  in  the  t-plane,  then 
from  (26),  (21)  and  (27): 
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u  *  -  o  + 


nK  sin  B 


0  e 


/■i 


dZ  _ 1 

dt  t,  t  -  t ' 


dZ 

nK  sin  B  > 


-v  •  „  dZl 

nK  sin  B  » 


c  <  t  <  0  ; 


0  <  t  <  e  ; 


nK  sin  B 


0  e 


c  0 


/-/  |i|t.  1“|t-t'|4t' 


The  evaluation  of  (45)  involves  calculating  a  Cauchy  Principal  Value  integral  which 
is  performed  by  subtracting  out  the  singularity  to  give 


u  «  -  a  + 


nK  sin  B 


0  e 


c  O' 


.  nK  sin  B  |  dZ  (t-c)(t-e) 
+ - - - |_  to - - 


Although  the  integrand  in  (47)  is  integrable,  for  numerical  purposes  a  first-order  estimate 
of  the  integral  in  the  region  of  the  logarithmic  singularity  is  made  separately. 


*  .?  dW  dt 

~u —  *  =  dZ  “  (u  "  lv)  dZ  ’ 


and  as  the  flow  is  conical 


4X  -  (<f  -  y^y  -  “  KU(I  '  yty  -  *$$)  * 


from  Euler's  theorem  on  homogeneous  functions. 

The  pressure  coefficient,  according  to  slender-body  theory,  is  given  by 

p  ‘  p~  .  vm  2  2*x  «  *  0 

P  iou  U2  U  IT 


i  .  (.f  -  2,« .  .(♦?  *  ;>-) 

ic  'K'  k2u  k2u2 


and  can  be  evaluated  from  the  velocity  components  given  in  (49)  and  (50) . 
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The  coefficient  of  normal  force  on  one  flap,  ,  then  follows  by  integrating  the 

jump  in  pressure  coefficient,  AC^  ,  over  the  flap,  and  non-dimensional ising  with  the  wing 
area  with  flaps  undeflected  to  give 


C 


NF 


where  £  is  measured  along  the  flap.  Using  the  fact  that  pressure  is  constant  along 
conical  rays  and  writing  £  »  ?s  gives 


(52) 


From  (19)  the  components  of  C^,  in  the  x  and  z  directions  are,  taking  leading 
order  terms,  -  sin  8  nKC^  and  cos  ®  cjjp  respectively.  The  coefficient  of  normal 
force  on  the  undeflected  part  of  the  wing  is  then  Cl3  -  2  cos  8  .  To  obtain  the 

drag  it  is  necessary  to  resolve  these  components  of  force  in  the  free-stream  direction  so 
that  the  drag  coefficient  (based  on  the  area  of  the  wing  with  undeflected  flaps)  at  the 
attachment  incidence  is  given  by  (for  the  slender-body  approximation  of  small  a  ) 


(\ 


2  cos  6  cNF)<*a  ♦  2  cos  6  "  2nK  sin  6  C^.  , 


or 


CD 

a 

7" 


La  “a  .  CNF 

— Y  “jf  ~  2n  sin  8 


(53) 


A  measure  of  the  lifting  performance  of  the  wing  is  then  given  by  the  lift- 
dependent  drag  factor,  x  *  where 


*acd  _ 

c.2  C.2  "k3"  \K2/ 


(54) 


and  where  here,  A  is  the  aspect  ratio  of  the  wing  with  undeflected  flaps. 

It  is  of  course  possible  to  obtain  the  force  on  the  undeflected  part  of  the  wing  by 
integrating  the  pressure  jump  over  this  region  in  a  similar  way  to  that  described  above: 
the  main  differences  are  that,  in  (45)  the  integral  is  no  longer  of  the  Cauchy  Principal 
Value  type;  (46)  becomes  v  ■  0  ;  and  in  (47)  no  logarithmic  singularity  is  involved  in 
the  integrand.  This  was  done  for  a  few  cases,  firstly  to  obtain  a  complete  picture  of 
the  pressure  distribution  on  the  wing  and  secondly  to  provide  a  check  on  the  method, 
since  the  total  force  in  the  z-direction  should  equal  the  lift  as  calculated  in 
section  4.4. 
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5  RESULTS 

The  geometry  of  the  wing  imposes  bounds  on  the  physically  attainable  parameter  space 
(n,8).  For  ft  <  0.5  the  opposite  flaps  will  meet  before  S  ■  it  .  In  fact  the  maximum 
angle  of  deflection  possible  for  a  given  ft  is 

8  *  cos  '( — ^ \  0  <  ft  <  0.5, 

“*  \n  -  l  / 


■  n  0.5  ^  n  ^  1.0. 

The  boundary  of  the  parameter  space  is  shown  in  Fig  3  together  with  the  region 
covered  numerically .  Difficulties  arose  in  the  convergence  of  the  Newton  scheme  for  find¬ 
ing  the  parameters  b,  c,  e,  f  in  the  transformation  near  the  (ii,8)  boundary  for 

n  <  0.5  .  This  was  due  to  1  and  F  coalescing  on  the  boundary,  so  that,  in  the  course 

of  the  iteration,  e  became  greater  than  f  ,  violating  (6). 

The  values  of  the  incidence  parameter  at  the  attachment  condition,  a^/K  ,  are 
shown  in  Fig  4  over  the  range  of  values  of  8  and  n  for  which  solutions  have  been 
obtained.  As  would  be  expected,  the  attachment  incidence  increases  as  the  deflection 
increases,  up  to  quite  large  deflection  angles.  However,  if  n  >  0.5  ,  so  that  the  flaps 
can  be  deflected  right  back  under  the  wing,  the  attachment  incidence  eventually  falls 
again  as  the  deflection  angle  approaches  ir  .  In  the  limit,  the  wing  is  again  a  flat 
plate,  giving  a  mathematical  check,  despite  the  physical  unreality  of  the  flow  model. 

The  condition  that  the  flow  attaches  at  the  leading  edge  of  the  flap  now  expresses  that 

the  flap  leading  edge  is  an  attachment  line  on  the  lower  surface  of  the  wing.  The  non- 

dimensional  coordinate,  nA  ,  of  the  attachment  line  He  the  distance  of  the  attachment 
line  from  the  centre  line  divided  by  the  local  semi-span),  on  the  lower  surface  of  a 
slender  flat-plate  delta  wing  is  given  by 


see  eg  Ref  9. 

The  semi-span  of  the  new  wing,  consisting  of  the  undeflected  portion  of  the  original 
wing,  is,  of  course,  ns  .  The  flap  leading  edge  is  s(l  -rj)  inboard  of  the  new  leading 
edge  (the  old  hinge  line)  and  so 


ns  -  s(>  -  ft) 


ns 


2ft  -  1 
ft 


Moreover  K  in  the  equations  above  must  be  replaced  by  tan  y’  ■  Kn  .  Then 
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-g-  “  n/l  -  ■  /(I  -  n) (3n  -  1)  for  n  >0.5  (55) 

or 

n  0.6  0.7  0.8  0.9 

a  /K  0.5657  0.5745  0.5292  0.4123  . 

a 

These  values  of  ofl/K  at  8  *  it  »  taken  with  the  computed  results  in  Fig  4,  imply  a  very 
rapid  variation  as  6  tends  to  it  ,  which  is  a  plausible  behaviour. 

The  practically  relevant  values  of  6  are  less  than  ir/2  ,  for  which  the  presenta¬ 
tion  in  Fig  4  is  confused.  Fig  5  shows  the  variation  of  <*a/$K  with  8  for  fixed  values 
of  n  •  The  values  at  8*0  correspond  to  the  fully-linearized  theory  (Appendix).  From 
equation  (A-l)  we  find  that,  for  the  linear  theory,  aa/8K  is  a  maximum  for  n  -  0.652  , 
taking  the  value  0.3572.  Fig  5  shows  that  for  8  >  0  ,  a  larger  maximum  is  obtained  at  a 
smaller  value  of  n  .  For  n  <0.8  ,  the  attachme  it  incidence  grows  more  rapidly  with  8 
than  the  linear  theory  predicts;  for  n  *  0.9  it  grows  less  rapidly.  If  n  ■  0  ,  so 
that  only  the  flaps  remain,  we  have  a  'caret  wing'  configuration.  As  n  0,  e  -*•  f  and 
in  (30)  two  singularities  of  the  integrand  coalesce,  the  integral  remaining  finite  for 
8  <  (tt/2)  ;  the  attachment  incidence  thus  tends  to  zero. 

The  comparison  with  linear  theory  is  brought  out  more  clearly  in  Fig  6,  where  the 
ratio  of  the  'exact'  attachment  incidence  to  its  linear  estimate  is  plotted  against  8 
for  fixed  values  of  n  .  For  inboard  positions  of  the  hinge  line,  large  errors  appear, 
but  for  the  most  practical  positions  at  80%  and  90%  of  the  semi-span,  the  linear  theory 
prediction  is  within  5%  for  deflection  angles  up  to  90°.  Even  for  a  hinge  line  at  60% 
semi-span,  which  is  probably  the  furthest  inboard  that  could  be  contemplated  on  a  real 
configuration,  the  error  is  within  10%  up  to  a  deflection  of  50°  and  only  exceeds  30%  at 
deflections  above  85°.  At  first  sight,  the  accuracy  achieved  by  the  linear  theory  is 
surprising.  However,  it  must  be  remembered  that  the  linear  result  is  obtained  by  applying 
the  boundary  conditions  on  that  part  of  the  plane  occupied  by  the  wing  with  the  flap 
undeflected.  This  is  not  the  procedure  adopted  in  the  application  of  linear  theory  to 
cambered  and  twisted  wings,  in  which  the  boundary  condition  is  applied  on  the  projection 
of  the  wing  onto  its  mean  plane.  If  the  conventional  procedure  were  applied,  it  would 
clearly  produce  a  substantial  underestimate  of  the  attachment  incidence  for  the  larger 
deflection  angles. 

2 

The  lift  coefficient  at  the  attachment  condition,  in  its  similarity  form  Cl^/K  , 
is  shown  in  Fig  7,  again  as  a  function  of  8  for  fixed  values  of  n  •  For  the  larger 
values  of  8  there  is  a  general  resemblance  to  the  curves  in  Fig  4.  In  particular  the 
same  rapid  variation  near  6  ■  it  arises.  The  values  for  6  *  it  follow  from  the 
familiar  expression  for  the  lift  of  a  flat  slender  wing. 


L 


2  2 

irpU  s  a  . 
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By  replacing  a  by  ns  and  introducing  S 


s2/K  and  a  from  (55) 
a 


n  0.6  0.7  0.8  0.9 

CLa/K2  1.280  1.769  2.128  2.098  . 


For  the  smaller  values  of  6  ,  the  lift  varies  more  rapidly  with  n  than  does  the  attach¬ 
ment  incidence. 

The  variation  for  the  smaller  values  of  0  is  brought  out  more  clearly  in  Fig  8, 

where  C,  /8K2  is  shown  for  0  <  8  <  1.5  .  The  values  at  0-0  correspond  to  the  fully- 
-a 

linearized  theory  given  by  equation  (A-3) .  It  is  easy  to  see  from  this  equation  that,  for 
a  given  deflection  angle,  the  linear  theory  predicts  a  maximum  value  of  C  at 
n  -  fay  -  0.816  .  The  present  calculations  show  that  a  hinge  position  near  802  of  the 
semi-span  is  close  to  the  optimum  for  deflection  angles  up  to  90  ,  though  little  is  lost 
by  any  choice  between  702  and  902. 

It  is  apparent  from  Fig  8  that  the  linearized  approximation,  in  which  the  curves  are 
replaced  by  horizontal  lines,  is  a  very  good  one,  particularly  for  the  three  most  outboard 
positions  of  the  hinge  line,  which  are  also  those  of  most  practical  interest.  The  points 
indicated  on  the  curves  show  the  extent  of  Che  range  for  which  the  linearized  result  is 
within  52  of  the  exact  value. 

The  lift-dependent  drag  factor,  x  »  i8  shown  in  Fig  9  for  the  three  most  outboard 
hinge  positions,  ie  n  -  0.7,  0.8  and  0.9,  and  for  flap  deflections  up  to  120°.  It  may 
be  noted  here  that  calculation  of  the  drag  caused  certain  problems,  not  least  being  that 
a  very  fine  subdivision  of  the  flap  was  necessary  in  the  numerical  evaluation  of  the 
integral  of  the  pressure  jump  across  the  flap  (equation  (52))  to  achieve  a  relatively 
low  accuracy  in  the  normal  force:  the  problem  worsened  as  deflection  and  flap  length 
increased.  This  difficulty  ultimately  limited  the  parameter  space  that  was  covered,  for 
the  lift-dependent  drag  factor,  because  the  computational  times  increased  beyond  accept¬ 
able  limits.  However,  the  most  physically  realistic  cases  fall  within  the  region  of  the 
(n,8)  space  covered  by  the  calculations. 

Again  comparison  is  drawn  between  the  linearized  boundary  condition  results  (dashed 
line)  given  by  equation  (A-8)  and  the  'exact'  boundary  conditions  (continuous  line),  with 
the  points  indicated  having  the  same  meaning  as  stated  above  for  the  lift.  Here,  also, 
the  linearized  approximation  is  seen  to  be  very  good  for  these  cases,  with  any  flap 
deflection  up  to  45°  being  within  the  52  error  band  for  all  three  hinge  positions.  Moving 
the  hinge  line  outboard  increases  the  range  for  which  the  error  in  the  approximation  is 
small. 

It  must  be  emphasized  that  all  coefficients  have  been  referred  to  the  area  of  the 
basic  wing  with  undeflected  flap,  and  it  is  the  aspect  ratio  of  this  wing  which  appears  in 


the  definition  of  x  •  Consequently,  the  solid  lines  in  Fig  9  show  the  reduction  in  lift¬ 
ing  efficiency  which  results  from  increasing  the  flap  deflection  on  a  given  wing  so  as  to 

maintain  the  attachment  condition  at  increasing  values  of  the  incidence  or  lift  coeffi¬ 
cient:  the  broken  lines  are  calculated  on  the  same  basis.  A  very  different  picture  arises 
if  the  flap  deflection  is  regarded  as  a  form  of  conical  camber  applied  to  a  wing  of  a  given 
planform,  ie  given  projected  area  in  the  plane  z  •*  0  .  If  x  is  written  in  terms  of  the 

drag,  D  ,  the  lift,  L  ,  and  the  semi-span,  s  ,  the  result  is 


2  2 

2irpU  s  D 

~u~ 


Consequently,  if  the  wing  remains  the  same,  with  the  same  forces  acting  on  it,  but  the 
span  is  defined  differently,  x  change*.  In  particular,  if  the  span  of  the  planform  is 

used  in  place  of  the  span  of  the  basic  wing  with  the  undeflected  flap,  x  is  reduced  by 

_  2  _ 
a  factor  (n  +  (1  -n)  cos  6)  •  For  the  rather  extreme  example  of  n  *  0.7,  B  ■  (it/2), 

this  means  reducing  the  value  of  1.71  shown  in  Fig  9  to  0.84.  This  value  is  comparable 

with  those  found  for  large  amounts  of  conical  camber  of  circular  arc  form  in  Ref  6. 

Figs  10,  II  and  12  show  the  pressure  distribution  over  the  whole  wing  (drawn  with 
the  flap  undeflected  for  ease  of  presentation)  for  the  cases  n  ■  0.7,  8  »  0.5  and 
n  3  0.8,  6  *  0.5  and  0  *>  1 . 1  .  As  stated  earlier,  the  suction  at  the  hinge-line  on 
the  upper  surface  is  infinite  although  no  localised  suction  force  is  present  for  flap 
deflections  below  0  »  tt  .  For  the  cases  with  the  same  flap  deflection  the  lifts  are 
similar  and  the  loading  can  be  seen  to  be  of  a  similar  magnitude.  For  the  larger  flap 
deflection  the  loading  both  on  the  undeflected  part  of  the  wing  and  on  the  flap  itself  is 
appreciably  larger  than  for  the  smaller  deflection  for  the  same  hinge-line  position. 

All  the  results  obtained  are  affected  by  the  limited  applicability  of  slender-body 
theory.  For  Mach  numbers  greater  than  unity,  the  linearized  theory  of  supersonic  flow 
provides  useful  guidance.  Equation  (A-2)  shows  that  the  lift  of  a  delta  wing  with  a 
slightly  deflected  leading-edge  flap  at  its  attachment  condition  increases  as  the  Mach 
number  increases  above  unity.  The  close  agreement  between  the  present  calculations  and 
the  fully-linearized  expression  (A-3)  suggests  that  (A-2)  is  likely  also  to  be  a  useful 
approximation  for  large  flap  deflections,  provided  the  leading  edges  are  not  too  close  to 
the  Mach  cone  from  the  apex.  For  subsonic  speeds  no  such  simple  result  is  possible;  since 
the  flow  over  the  conical  configuration  will  no  longer  be  conical,  there  will  be  no 
incidence  at  which  the  flow  will  be  attached  along  the  whole  leading  edge.  Practically 
this  seems  to  be  unimportant.  For  instance,  a  wing  of  aspect  ratio  0.75,  with  a  gothic 
planform,  8.25)5  thick  on  the  centre  line,  with  camber  and  twist  designed  by  slender  theory 
for  attached  flow  at  a  lift  coefficient  of  0.1,  was  found*^  to  show  no  signs  of  leading- 
edge  separation  at  low  speeds  at  lift  coefficients  of  0.09  and  0.15. 

6  CONCLUSIONS 

The  application  of  slender-body  theory  to  a  delta  wing  furnished  with  a  conical 
leading-edge  flap,  deflected  so  as  to  maintain  attached  flow  at  the  leading  edge,  shows 
the  following: 


18 


(1)  Increasing  the  angle  of  deflection,  measured  in  the  cross-flow  plane,  up  to  90°  and 
beyond  produces  an  increasing  angle  of  incidence  for  attachment  and  an  increasing  lift 
coefficient  at  the  attachment  condition.  At  the  same  time,  there  is  a  decrease  in  the 
lifting  efficiency,  as  measured  by  a  lift-dependent  drag  factor  based  on  the  planform  of 
the  wing  with  the  flap  undeflected. 

(2)  The  fully-linearized  slender  thin-wing  theory  agrees  fairly  well  with  the  slender- 
body  theory  for  the  practical  cases  of  outboard  hinge-lines.  The  lift  at  attachment  is 
predicted  particularly  closely  by  the  linear  theory  for  hinge  lines  at  70Z,  80Z  and  90Z  of 
the  semi-span,  where  the  error  is  less  than  SZ  for  deflection  angles  up  to  75°.  The  lift- 
dependent  drag  factor  is  within  the  same  tolerance  for  deflection  angles  up  to  45°,  for 
the  saaie  hinge-line  positions. 

(3)  The  largest  lift  coefficients  at  attachment  for  a  given  deflection  angle  are 

produced  with  the  hinge-line  near  80Z  of  the  semi-span,  but  the  degradation  at  70Z  and  90Z 

is  small.  For  this  hinge-line  position,  Che  lift  coefficient  at  attachment  is  about 
2 

1.5  BK  ,  where  K  is  the  tangent  of  the  semi -apex  angle  and  0  is  the  deflection  angle 
(in  radians),  for  8  <  (tt/2)  . 
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Appendix 

LINEARIZED  BOUNDARY  CONDITIONS 

Smith  and  Mangier11  calculated  the  supersonic  flow  over  a  delta  wing  with  hinged 
flaps  at  the  subsonic  leading  edges  using  linearized  boundary  conditions.  Their  result 
(in  closed  form)  for  the  attachment  incidence  is 


a 


a 


n 


J  U- 

rl  -  a  n 


where  here  a 
number,  and 


tan  y/tan  p 


where 


p  ■  cosec  M 


is  the  Mach  angle. 


M  being  the  Mach 


n 


n^a 


2-1 
a  n 


is  the  complete  elliptical  integral  of  the  third  kind.  In  this  case  £  -  -nK8  .  The 
slender  wing  corresponds  to  the  limit  a  -  0  .  Taking  this  limit  we  find  that 


a  ~  I  — 

a  _  2n  cos  n 

Kg  it 


The  lift  coefficient  is  stated  to  be 


which  yields,  on  taking  the  same  limit 


(A-l) 


(A-2) 


(A-3) 


We  now  show  that  these  expressions,  (A-l)  and  (A-3),  are  obtained  from  the  present 
treatment  by  retaining  terms  of  leading  order  in  6  only.  Equation  (30)  shows  that 
aa/K  is  proportional  to  8  when  6  is  small,  so  the  integral  in  (30)  may  be  evaluated 
for  8  ■  0  if  only  the  leading  term  is  needed.  By  the  same  argument,  in  (44)  it  is 
adequate  to  evaluate  the  integral  and  the  coefficient  y_(  for  8  ■  0  .  Then  by  (35) 


y_l  *  -  |  (b  -  f)2 


(A-4) 


where  the  overbar  denotes  the  quantity  evaluated  for  6  ■  0  .  To  find  b,  c,  e,  f  we 
put  8  •  0  in  (14)  to  (17).  Adding  (16)  and  (17)  and  subtracting  (14)  and  (15)  gives 
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With  t  « 


and  so  b 


With  t  - 


from  which 


The 


and  so 


to  leading 


/ 


t _ 

b) (f  -  t) 


dt 


0  . 


-2-2 
b  cos  6  +  f  sin  9 


the  integral  reduces  to 


it/ 2 

2  J  (b  cos20  +  f  sin20)d9  -  j  (b  +  f) 
0 


-  -f  .  Now  adding  (16)  and  (17),  again  with  6  "  0  ,  gives 

f 


f  sin  0  ,  the  integral  reduces  to  f  ,  so  f  «  1  ,  Equation  (17)  reduces  to 


/ 


n 


By  inspection  of  (15)  and  (16) 


c  -  -  e  . 


integral  in  (30)  becomes 


.  -1- 
sin  e  - 


.  -I_ 
sin  c 


-1_ 
cos  n 


aa  2g?i  cos  1  n 

K  7T 


(A-5) 


order.  The  integral  in  (44)  becomes 


/ 


dt 


-1_ 

cos  n 


and  so,  using  (A-3)  we  have 
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LIST  OF  SYMBOLS 

constant  in  transformation  also  aspect  ratio  of  ving  with  undeflected  flap 

-  tan  v/tan  y  (see  Appendix) 

constants  in  expansion  of  W  in  terms  of  Z 
coordinates  of  B,  C,  E,  F 

lift  coefficient  based  on  total  area  of  wing  and  flaps 
lift  coefficient  at  incidence  for  attachment,  a# 

normal  force  coefficient  on  flap 
pressure  coefficient 

drag  coefficient  at  incidence  for  attachment 

constants  in  expansion  of  t  in  terms  of  Z 
drag 

-  /n 

”  tan  y 

lift 

lift  at  attachment  incidence,  a 

’  a 

Mach  number  (see  Appendix) 
source  strength 

distance  and  non-dimensional  distance  normal  to  contour 

unit  outward  normal  to  wing 

wing  area  with  flaps  undeflected 

semi-span  of  cross-section  with  flaps  undeflected 

complex  coordinate  in  transformed  plane 

free  stream  velocity 

complex  velocity  in  t-plane  (non-dimensional) 

complex  potentials 

cartesian  coordinates 

••  y/s,  z/s,  non-dimensional  coordinates 

-  (y  +  iz)/s,  non-dimensional  complex  coordinate  in  cross-flow  plane 
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LIST  OP  SYMBOLS  (concluded) 


incidence 

incidence  for  attachaent 

flap  deflection  in  cross-flow  plane 

semi-apex  angle  of  wing  with  flap  undeflected 

semi-apex  angle  between  hinge  lines 

constants  in  expansion  of  Z  in  terms  of  t 

non-dimensional  coordinate  of  attachment  line,  equation  (SS) 

conical  coordinate  of  the  flap  hinge-line 

»  cosec  *M  ,  the  Mach  angle  (see  Appendix) 

density 

velocity  potential 
cross-flow  potential 
-  <{>/Us  ,  non-dimensional  potential 
stream  function  corresponding  to  $ 
lift-dependent  drag  factor,  equation  (54) 
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